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ABSTRACT 

We study the role of major and minor mergers in the mass growth of lu- 
minous red galaxies. We present small-scale (0.01 < r < 8 h^ 1 Mpc) pro- 
jected cross-correlation functions of 23043 luminous early-type galaxies from 
the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy (LRG) sample 
(0.16 < z < 0.30, [Mo.3j — 5 \og 1Q h] « —22.75 mag) with all their companions 
in the SDSS imaging sample, split into color and luminosity subsamples with 
[Mo.3j — 5 log 10 h] < —18 mag. We de-project the two-dimensional functions to 
obtain three-dimensional real-space LRG-galaxy cross-correlation functions for 
each companion subsample. We find that the cross-correlation functions are not 
purely power-law and that there is a clear "one-halo" to "two-halo" transition 
near 1 h~ l Mpc. We convert these results into close pair statistics and estimate 
the LRG accretion rate from each companion galaxy subsample using timescales 
from dynamical friction arguments for each subsample of the companions. We 
find that the accretion onto LRGs is dominated by dry mergers of galaxies more 
luminous than L*. We integrate the luminosity accretion rate from mergers over 
all companion galaxy subsamples and find that LRGs are growing by [1.7 ± 0.1] 
percent per Gyr, on average, from merger activity at redshift z ~ 0.25. This rate 
is almost certainly an over-estimate because we have assumed that all close pairs 
are merging as quickly as dynamical friction allows; nonetheless it is on the low 
side of the panoply of measurements in the literature, and lower than any rate 
predicted from theory. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution - 
galaxies: interactions — large-scale structure of universe — methods: statistical 
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Introduction 



In the current paradigm for galaxy formation, the massive dark-matter halos in which 
galaxies reside have assembled throughout cosmic time by accretion and merging of smaller 
parts. There remains, however, substantial uncertainty about the formation and evolu- 
tion of the galaxies in those halos, and the evolution of the stars and gas that compose 
them. A particularly important galaxy subpopulation in this research context is the lumi- 
nous end of the red sequence of galaxies — consisting of concentrated, smooth, and typically 



ellipt i cal galaxies that preferentially live in the densest regions of the Universe JSandag' 



1972; 



Schneider et al.lll983l : iRoberts fc Haynedll994j ; iPostman fc Lauerlll995l ; iBlanton et al. 



2003 bJ). Such galaxies account for a large fraction of the to tal stellar mass in the Universe 



( IHogg et al.l l2002bl ; iRudnick et al.ll2006l ; iBrown et al.l 120071 ). and thus understanding their 
formation is critical to understanding galaxy formation in general. This red sequence is 
clearly separated from the blue sequence of galaxies, which are typically lower in mass, more 
star-forming, gas-rich, morphologically "sp i ral," and preferent i ally populate isolat ed regions 
( Istrateva et aDboOll ; IBlanton et aDl2003bl : iBaldrv et aPl2004j ; iBalogh et aDl2004h . 



With uniform spectra, deep absorption lines, highly clustered distribution ( "bias" around 
2), and high luminosities (absolute magnitudes around —23 mag), th e Luminous Red Galax - 
ies (LRGs) are excellent tracers of the density field on large scales (jEisenstein et al.ll200ll ). 
For this reason, surveys of LRGs wer e the first to conc lusively demonstrate the homogene- 
ity of the Universe at low redshifts (IHogg et al.l 120051) and to dete ct the baryon acoustic 
oscillation feature in the correlation function (jEisenstein et al.l 120051 ). Their continuing im- 
portance to understanding fundamental cosmology underlines the need to better understand 
their nature. 

Because luminous red galaxies are typically supported by velocity dispersion and not 
orderly rotation, it has been h ypothesized that they form from mergers of two or more 



smaller galaxies (iToomrd 119771 ). Such mergers have been shown in numerical simulations 



cals 


Negroponte & White 


1983; 


Barnes & Hernauist 


1996; 


Naab & Burkert 


2003; 


Cox et al. 


2006 


). In addition to such "major" mergers, LRGs could in principle grow over time from 



an accumulation of smaller mergers. I f at least one of the merging galaxies is gas-rich, it 
often shows a large star- formation rate (IBarton et al.ll2000l ; lLambas et al.ll2003l ; iNikolic et al. 



2004 ; ISmith et al.l 120071 ) . and indeed a significant fraction of the total star-formation in the 
Universe may occur in such events. However, there is also a population of mergers of two red 
galaxies, in which no star-formation occurs, an event usually referred to as a "dry merger" 
(feell et al.l [2006b! : Ivan Dokkumll2005h . 



A growing consensus of groups studying the high redshift Universe find that the luminous 
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red galaxies appear to grow in stellar mass over time 


(Bell et al. 


2004; 


Willmer et al. 


2005; 


Blanton 


2006 




Wake et al. 


2006: 


Brown et al. 


2007 


: ] 


7 aber et al. 


2007 


). This growth can 



occur in several ways. First, in principle they may have ongoing star-formation. This is very 
unlikely, since the red galaxies show little signs of such star-formation. Second, luminous 
blue galaxies may transform to red galaxies. This is also unlikely at the very luminous end, 
since the number density of luminous blue galaxies is far lower than that of red galaxies, 
even at high redshift. Third, and most likely, the luminous red galaxies may grow through 
either major or minor mergers. 

If the change in the galaxy population over time is to be explained at least partly 
by mergers, we must be able to find these mergers before or as they occur in appropriate 
numbers. "Instantaneous" studies of the merger rate are complementary to, and must be 
consistent with, the global studies in the change of the galaxy population. To evaluate 
whether mergers can explain these changes, the work presented here is designed to measure 
the accretion or merger rate of companion galaxies into LRGs at redshift z ~ 0.25. 

It is also the case that in the CDM paradigm for structure formation, galaxies reside 
in mass concentrations that are built from merging and accretion of smaller concentrations 
over cosmic time. It is unavoidable that this mergi ng in the dark sector is associated, at 
some level, with m erging of observable galaxies (e.g. 
Conrov et aDboOTt ). 



Murali et al. 2002; Mailer et al. 2006 



There are many galaxy-galaxy merg er rate estimates in t he literature, which involve 



ident i fying a class of pre-merger close pairs (ICarlberg et al. 



199 



9; 



2006a 



Carlberg et al. 



2000 



De Propris et al. 



Patton et al. 2000 



2007) 



Lin et al 



1994 



Patton et a. 



20041 Masiedi et al 



. 1997 



2006 



van Dok 



cum et al. 



Bell et al. 



a class of post-merger galaxies based on star-formation in- 
dicators (jQuintero et al.ll2004f). or a class of cur r ently merging sources based on disturbed 
or merging morphologies jAbraham et aL 1996 ; Conselice et al. 2003 ; van Dokkum 2005 



Lotz et al. 



20061 : be Propris et al.ll2007h . In each case, the estimate of the merger rate pro- 



ceeds by estimating the abundance of the class, some time interval over which they remain 
identifiably part of that class. The measurements of the LRG accretion rate presented here 
also follow this methodology. However, our measurements are more reliable than most previ- 
ous measurements for a number of reasons. The first is that the LRGs form a very uniform, 
very massive population, as described above, and therefore dynamical times and dynamical- 
friction times relevant to close pairs are straightforward to estimate. The second is that we 
make maximal assumptions so as to put a strict upper limit on the accretion rate. This is 
interesting, because the upper limit we determine is on the low side of existing predictions 
and measurements. The third is that we use a technique for measuring the mean number 
of close pairs in real space, with no contamination by projected pairs, so our pre-merger 
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candidate list is clean of such interlopers (in a statistical sense). 

Our technique for measuring the close pairs builds on previous work (IMasjedi et al. hood ) 



in which we showed that we could measure projected correlation functions on extremely small 
scales (kpc to Mpc scales) without the need for complete spectroscopic samples. We used this 
method to overcome the fiber-collision incompleteness of the LRG sample in the SDSS. In 
addition, we showed that the clustering signal so measured can be deprojected and integrated 
to deduce close pair statistics and therefore a rate of merger events among the galaxies in 
the sample. We found that LRG LRG mergers are extremely rare events and do not play a 
significant role in the growth of these galaxies, at least at low redshifts. 

In this paper we expand this technique to measure not only the auto-correlation func- 
tion of a set of galaxies, but the cross-correlation function of two different galaxy sets, only 
one of which requires spectroscopic information. We choose LRGs as our primary spectro- 
scopic sample and we cross-correlate them with distinct subsamples of companion galaxies 
over a range of luminosities and colors. We convert these results into an accretion rate of 
luminosity — within each subsample — into LRGs. We can estimate a total accretion rate by 
this method with unprecedented precision. 

Throughout this paper, all distances are comoving, calculated for a cosmological world 
model with (f2 m , = (0.3, 0.7) and Hubble constant parameterized by H Q = 100 h km s _1 Mpc 
All magnitudes are AB. 



-i 



Data 



The SDSS (IStoughton et al.ll2002l ; lAbazajian et al.ll2003l ; lAbazaiian et al.ll2004f) has per- 



formed an imaging and spectroscop ic survey of ~ 10 4 square degrees (IFukugita et a 



image processi ng (ILupton et al 



ric ca libration (jSmith et al 



Gunn et al.lll998t iGunn et al.l 120061). Automated, real-time monitoring 



2001 



2002 



Stoug 



Ivezic et al 



iton et al.l 
200i 



2002 



Tucker et al 



Pier et al. 
20061 : 



Hogg 



et al 



1996 



20011 ) . 



2003 ), photomet- 



Padmanabhan et al. 



2007 ). galaxy target selection for spectroscopy (Strauss et al. 20021 ; lEisenstein et al. 2001 ) 
design of spectroscopic plates (IBlanton et al.ll2003al ). and spectroscopic reductions have pro- 
duced enormous, very uniform samples. Of the vario us SDSS subsamples, the one that 



uniformly maps the largest volume is the LRG sample (lEisenstein et al.ll200ll ). 
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2.1. Spectroscopic subsample 



The spectroscopic LRG sample is constructed from color-magnitude cuts in g, r, and i 
bands to select galaxies that are likely to be luminous early-type galaxies at redshifts between 
0.15 and 0.5. The selection is highly efficient and the redshift success rate is excellent. The 
sample is constructed to be close to volume-limited up to z = 0.36, with a dropoff in density 
toward z = 0.5. Because the LRG sample is so uniform, and because it oc cupies such a larg e 
volume, we have used it to demonstrate the homoge neity of the Universe (|Hogg et al.ll2005l ). 
to locate the baryon acoustic feature at low redshift (Eisenstein et al.ll2005l). and to measure 



clustering at intermediate and small scales (IZehavi et al 



2005 



Masiedi et all 120061 ) . 



This study uses a spectroscopic sample drawn from NYU LSS sample 14 (IBlanton et al. 



20051 ). This covers 3,836 square degrees and contains 55,000 LRGs with redshifts 0.16 < 
z < 0.47. The subsample of LRGs used in this paper has luminosity and redshift ranges 
of -23.2 < [Mo.3 g - 5 log 10 /i] < -21.2 mag and 0.16 < z < 0.30. We restrict our sam- 
ple to z < 0.30 to allow measurement of cross-correlations between LRGs and much less 
lumin ous companions. Th e LRG absolute magnitudes include Gal actic extinction correc- 
tions (ISchlegel et al.lll998l ). K corrections (IBlanton fc Roweid 120071 ) and passive evolution 
corrections (the latter were applied only for the purposes of selecting a sample that does 
not substantially change with redshift). These cuts left 23043 LRGs in our spectroscopic 
subsample. 



pectroscopic fi 


3er pluj 


(IBlanton et al. 


2003a) 



lar completeness is characterized for each unique region of overlapping spectroscopic plates 
("sector") on the sky. An operational constraint of SDSS spectrographs, however, is that 
the physical size of the fiber coupling forces the angular separation of targets on any indi- 
vidual spectroscopic plate to be larger than 55 arcsec. This "fiber collision" constraint is 
partly reduced by having roughly 40 percent of the sky covered by overlapping plates, but it 
still results in ~ 7 percent of targeted galaxies not having measured redshifts. Because this 
project involves cross-correlating spetroscopic and imaging objects, this fiber collision limit 
only comes into our analysis in our weighting scheme to account for incompleteness; it does 
not affect our pair counts directly. 

For each galaxy j in the spectroscopic subsample we compute a weight Pj that accounts 
statistically for the spectroscopic incompleteness coming from fiber collisions. We calculate 
this weight by running a two-dimensional friends-of-friends grouping algorithm on the SDSS 
target parent sample in sam plel4, with a 55 arcs ec linking length. This procedure emulates 
the SDSS tiling algorithm (IBlanton et al.ll2003al ). Within each "collision group" made by 
the friends-of-friends algorithm, we find the number of objects with measured spectroscopic 
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redshifts and divide by the total number. The inverse of this ratio is a weight pj assigned to 
each spectroscopic LRG to account for survey incompleteness. 

We have created large catalogs of randomly distributed points based on the SDSS sub- 
sample angular and radial (redshift distribution) models. These catalogs match the redshift 
distribution of the LRGs and are isotropic within the survey region. These catalogs allow us 
to check the survey completeness of any given volume and provide a homogeneous baseline 
(expected numbers) for the tests that follow. 

For each random point j we compute a weight fj that accounts for the incompleteness 
of the spectroscopic survey in that point's region of the sky not due to fiber collision but 
due to all the other selection effects in the survey. The sample 14 package provides the 
angular geometry of the spectroscopic survey expressed in terms of spherical polygons. The 
geometry is complicated: the spectroscopic plates are circular and overlap, while the imaging 
is in long strips on the sky, and there are overlap regions for some plates that have not yet 
been observed. The resulting spherical polygons track all these effects and characterize the 
geometry in terms of "sectors", each being a unique region of overlapping spectroscopic 
plates. In each sector, we count the number of possible targets (LRG, Main, and Quasar), 
excluding those missed because of fiber collisions, and the number of these whose redshifts 
were determined. We weight the points in the random catalog matched to the spectroscopic 
LRGs by the inverse of the ratio of these numbers (fj). In truth, the priority of all targets 
are not equal, such that LRGs always "lose" to quasar candidates, but the LRG priority 
is equal to that of the dominant MAIN targets. Only about 12 percent of the fibers are 
assigned to quasars, hence quasar-LRG collisions are rare and this priority bias is small. 

We are required to treat the fiber collision incompleteness pj and overall incompleteness 
fj factors separately, because the former is strongly correlated with LRG environment, and 
there can be physical differences between LRGs in high and low density environments. 



2.2. Imaging subsamples 



For our imaging data we use the full imaging sample of the SDSS imaging catalog in the 
DR4plus footprint, which completely covers samplel4 and is equivalent to the SDSS Data 
Release 5. After correcting for Galactic extinction, we applied an i-band apparent magnitude 
cut of rrii < 21 mag. This cut guarantees completeness in the redshift range z < 0.30, K- 
corrected absolute magnitude range [Mo.3, — 5 log 10 h] < —18 mag. The °' 3 g, °' 3 r, and °' 3 i 
bandpasses are the SDSS g, r, and i bandpasses shifted bl ueward by a facto r of 1.3 so that 



K corrections for galaxies at z — 0.3 become trivial (e.g., iHogg et al.ll2002al ; iBlanton et al. 
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2003bl ). We applied an additional surface-brightness cut of fi < 28 mag in 1 arcsec 2 in g, r 
and i. The surface-brightness cut is far below the SDSS detection limit; it cleans the data 
of the obvious mis-measurements of the Petrosian aperture and other extended-source data 
artifacts. 

We cannot .fT-correct individual galaxies in the imaging subsample once and for all, 
because we do not have spectroscopic redshifts for them, but each time we consider a pair 
of galaxies, one from the spectroscopic subsample and one from the imaging subsample, we 
fictitiously assign the spectroscopic redshift to the imaging galaxy. This allows us to calculate 
for each imaging galaxy in each spectroscopic-imaging pair a "temporary" K-corrected °' 3 i- 
band absolute magnitude and [ 0,3 g — 0,3 r] color for the purposes of that pair. We discard 
these values and compute new ones when the imaging galaxy is used in another pair with 
another spectroscopic galaxy. 



We calculated the K-corrections using the code kcorrect (IBlanton fc Roweid 120071 ) . 
This code is accurate but too slow to calculate the if-corrections individually for the number 
of pairs (~ 10 9 ) found in the cross-correlations. To save time, we computed the K-correction 
on a grid of colors in advance. We took galaxies from the SDSS Main Sample as representative 
of all galaxy types. We computed their K-corrections on a grid of redshifts between 0.16 and 
0.30 (the redshift limits of our spectroscopic subsample). We saved the mean i^-correction 
in a grid of observed [g — r] color, [r — i] color, and redshift. Thereafter we interpolated 
this cube when calculating the i^-correction for an galaxy in any imaging subsample. This 
speeds up the i^-correction procedure immensely and only introduces percent-level errors in 
the results. 

We have created large catalogs of randomly distributed points, with the angular distri- 
bution of the imaging data subsamples. 



3. Method and results 

As we describe in this section, we cross-correlate galaxies in a spectroscopic subsample s 
(of LRGs in this case) with galaxies in an imaging subsample % to obtain the real-space, pro- 
jected cross-correlation function w si (r p ) as a function of tangential projected separation r p . 
We de-project this projected cross-correlation function to obtain the true, three-dimensional, 
real-space cross-correlation function £ S j(r) as a function of real-space separation r. We use 
this three-dimensional cross-correlation function and dynamical arguments to place limits on 
the accretion rate of objects from subsample i into objects from subsample s, and therefore 
the mass growth rate of LRGs. 
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3.1. Projected cross-correlation function 

In analogy to the definition of auto-correlation function, the three-dimensional real- 
space cross- correlation function £ S j(r) of two subsamples of galaxies s and i, is defined as the 
excess probability of finding a galaxy from subsample s at a distance r from a galaxy from 
the subsample i, relative to the "null" Poisson prediction. If we take two small comoving 
volumes dV s and dVi, in which we look for galaxies from subsamples s and i respectively, 
separated by a distance r, the expected number of pairs dN S i with one galaxy coming from 
subsample s and the other from subsample i is: 

dN ai = n a m[l + ^ ai (r)]dV a dV i , (1) 

where n s and rii are the three-dimensional comoving number densities of subsamples s and 
i respectively 

The projected two-dimensional cross-correlation function w S i(r p ) is related to the three- 
dimensional real-space correlation function £«i(r) by a projection over the component tt of 
the separation along the line of sight 



w si (r p ) = J d7r6i(yr* + 7T 2 J . (2) 

The two-dimensional function w S i(r p ) has dimensions of length. Because in practice the 
correlation function £ S j(r) is very large at small scales, the integral is dominated by scales 
7i < r p . Observationally, w S i(r p ) is much more accessible than £ s j(r), because the radial 
component of the separation is never well measured (and not measured at all in the work 
presented here). 



Following the approach we have used previously (jMasjedi et al.l 120061 ) , we measure 
w S i(r p ) schematically as a difference of two ratios: 

D s Di R s Di 

n iWsi {r p ) = ——-—— , (3) 
U s tii rC s Ki 

where n» is the average comoving three-dimensional volume density of the imaging subsample, 
the symbols D s and Di represent the spectroscopic and imaging data subsamples, and R s 
and Ri represent the random catalogs matched to the spectroscopic and imaging subsamples 
respectively. The product of a volume density and a length, niW si (r p ) has dimensions of 
inverse (comoving) area. In a rough sense, the first term on the right-hand side of equation 
([3]) measures the abundance of pairs, and the second term subtracts the mean background 
level. The procedure described here has been tes ted with simulation s and shown to deliver 



an unbiased measure of the correlation function (jMasjedi et al.ll2006f ). 
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The main difference between equation ([3]) and our previous work (IMasjedi et al.l 120061 ) 
is that the number density n.j that enters on the left-hand side is the number density of 
the imaging subsample, which we cannot determine explicitly within this data set since, by 
construction, the imaging subsample has no (or few) spectroscopic redshifts. We can only 
measure robustly the product riiW S i{r p ). Fortunately, for the purposes of estimating the 
merger rate, we need only this product, and not either quantity separately. 

The right-hand side of Equation includes a spectroscopic-imaging pair-count factor 
D s D f 

52 p i 

D s D l = ^— , (4) 

keo s 

where the top sum is over pairs j with one member taken from the spectroscopic subsample 
and one from the imaging subsample in some bin of transverse radii r p , the bottom sum is over 
galaxies k from the spectroscopic subsample, and pj is the weight given to the spectroscopic 
galaxy in pair j that accounts for fiber-collision incompleteness as described above. Being a 
sum of dimensionless weights, this factor D s Di is dimensionless. 

There is a spectroscopic-random pair-count factor D s Ri\ 

52 p^> 

D,R, = ,„ , (5) 



52 P k 



dtt 

dA 



dN 

dn 



where the top sum is over pairs j with one member taken from the spectroscopic subsample 
and one taken from the random catalog matched to the imaging subsample, the bottom sum 
is over galaxies k in the spectroscopic subsample, [^] , is the inverse square of the comoving 
distance to the spectroscopic galaxy k, and 4?| is the number density of the random imaging 
catalog per solid angle. The product of these two derivatives gives the average number of 
random imaging objects per unit comoving area around each spectroscopic galaxy, so this 
factor D s Ri has dimensions of comoving area. 

There is a random-imaging pair-count factor RsDf. 



E n 



k£R s 



10 



where the top sum is over pairs j with one member taken from the random catalog matched 
to the spectroscopic subsample and one taken from the imaging subsample, the bottom sum 
is over points in the random catalog matched to the spectroscopic subsample, and fj is 
the weight given to the random point in pair j that accounts for the incompleteness of the 
spectroscopic survey in that point's region of the sky not due to fiber collisions, as described 
above. This factor R s Di is dimensionless. 

There is a random-random pair-count term R s Rf. 



similar to the above but for pairs j with one member taken from the random catalog matched 
to the spectroscopic subsample and one taken from the random catalog matched to the 
imaging subsample. This factor R s Ri has dimensions of comoving area. 

In measuring the four factors on the right-hand side of Equation ([3]), we have used 2 
imaging-galaxy color bins separating the blue and red imaging galaxies in [ a3 g — 03 r] color, 
20 imaging-galaxy luminosity bins, which we choose to cover the range —24 < [Mo.3j — 
5 log 10 h] < —18 mag but have roughly the same number of imaging galaxies in each, and 15 
bins in projected transverse separation r p between the LRG and the accompanying galaxy, 
covering the range of 0.01 < r p < 8 h^ 1 Mpc with logarithmic spacing. 

In practice, to compute the factors, we bin all the specroscopic-imaging pairs according 
to the imaging galaxy color, the imaging galaxy luminosity, and the comoving projected sep- 
aration r p of the pair. As described above, the K corrections and separations are computed 
for each pair using the redshift of the spectroscopic galaxy. We perform the sums given in 
Equations (j4]) through ([7]) in each bin separately and thereby construct the w S i estimator 
given in Equation ([3]). 

Figures [T] and [2] show the measurements of ^ u> S i(r p ) for red and blue companion galax- 
ies respectively. We have combined the 20 luminosity bins into 5 to simplify the figures. The 
error bars are estimated using jackknife resampling covariance matrix with 100 subsamples 
made contiguous and compact on the sky (based on SDSS "targetting chunks") to be as 
conservative as possible with regards to correlated calibration and selection errors. Note 
that the error bars for each subsample are smallest on kpc scales and become larger for both 
the smaller and larger scales. On smaller scales this is due to shot noise; the smaller the sep- 
arations the fewer the pair counts. On scales larger than a few h~ l kpc, the errors grow both 
due to cosmic variance and the fact that our method becomes more and more vulnerable to 



R s Ri 




(7) 



interlopers on larger scales where the clustering power is weaker and background subtraction 
is more noisy. These effects generate high correlations among the errors of different bins, 
and explains the smoothness of the curves in Figures [T] and [2] despite the large uncertainties 
in each bin. 



3.2. Three-dimensional statistics 

Under the assumption of spherical symmetry, the two-dimensional, projected cross- 
correlation function w S i(r p ) can be "deprojected" into the three-dimensional, real-space cor- 
relation function £ S j(r): 

n iU r) = -- I" , (8) 

71 Jr V ?"p - r 2 dr p 

where we have kept this deprojection in terms of the measureable product niW si {r p ). 

The correlation function can be converted to pair counts; similarly the cross-correlation 
function can be converted into the mean number of galaxies from a specific imaging 
subsample i within a given small three-dimensional separation r c i ose of a member of the 
spectroscopic subsample s: 

/"^close /"^closc 

Ni^A-KUi r 2 dr [1 + £«(/•)] ^4tt / r 2 dr [m £ si (r)] , (9) 

Jo Jo 

where we have used the fact that at small scales £ S i ^> 1 and the term in brackets in the 
approximate expression is the quantity we can de-project from the two-dimensional projected 
cross-correlation function niW S i(r p ). So we can measure the close pair fraction for every 
subsample for which we can measure the cross- correlation function. 



3.3. Merger rate 

Conversion of a pair fraction measurement into a merger rate requires a time-scale t m erge,i 
over which the mean galaxy from imaging subsample i within separation r c i ose will merge 
with the mean LRG from subsample s. The merger rate estimate Tj of galaxies from sample 
i into galaxies from sample s per galaxy (from s) per unit time is 
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and the mean fractional rate of growth of luminosity of a galaxies from subsample s from 
accretion of galaxies from subsample i is 



(L s ) 



dLs 
dt 



"merge,* (L s ) 



where (L s ) is the mean luminosity of galaxies from subsample s and (Lj) is the mean lumi- 
nosity of galaxies from subsample i. 

The shortest conceivable merger time t mer ge,; estimate (which produces the largest con- 
ceivable estimate of the merger rate) is the orbital time £ rbit- A more realistic estimate is 
a time td yn ,i based on dynamical friction. But in principal all of these times can be under- 
estimates (and hence any merger rate based on close pairs can be an overestimate) because 
there is undoubtedly a large number of close pairs that will not merge on any short timescale. 
In what follows, we present the orbital time t rbit and dynamical friction time td yQ ,i as two 
options, but then interpret our merger rate estimates as upper limits. 

All of these merger rate estimates depend, in principle, on the radius r c i osc inside of 
which we have counted close companions. However, over the range of interest in Figures [1] 
and El w S i(r p ) scales (something) like r p , £ s i(r) scales (something) like r 2 and N scales 
(something) like r c i ose . Similarly, both time-scales (orbital and dynamical-friction) scale like 
r ciose- For this reason, the inferred merger and accretion rates (which are based on ratios of 
Ni with the timescales) do not depend strongly on the choice of r c i ose . 

The average orbital velocity for a companion around a more massive galaxy with velocity 
dispersion a v is roughly 1.5 times the velocity dispersion, so 

torbit^^^ ■ (12) 

1.5 a v 

This is the shortest conceivable mean merger time (we have included the factor of 1.5 to be 
conservative). The fractional luminosity accretion rate estimate for this assumed time-scale 
is shown with dashed lines in Figure [3] as a function of the luminosity of imaging subsample 
i for red and blue imaging galaxies. The per-subsample merger rates have been divided by 
the absolute-magnitude bin width so that the total fractional accretion rate is the area under 
(integral of) the curves. 

The Chandrasekhar approximation to dynamical friction is longer than the dynamical 
time by a factor roughly equal to the ratio of the mass of the heavier galaxy to the lighter 
one. This approximation may a ctually be an underestimat e of the total merger time found 



in explicit iV-body simulations (jBoylan-Kolchin et al.l 120071 ). which serves to strengthen the 



interpretation of our merger rate estimate as an upper limit. For our case, the approximation 



13 



becomes 

(m s ) 

tdvnA ^orbit "7 T > (.1") 

(m { ) 

where (m s ) and (m,) are the averages of the masses of the spectroscopic and imaging sub- 
samples respectively, and we have assumed (m s ) > (mi). Keeping things observational, we 
do not try to measure masses for galaxies in this work. Instead, we make the naive assump- 
tion that a galaxy's mass is directly proportional to its 3 i-band luminosity and therefore 
we use the ratio of the luminosities instead of the masses. The 03 i-band luminosity is very 
close to the rest frame r-band luminosity. This assumption works fairly well for the mass 
ratios of red imaging galaxies to the spectroscopic galaxies (which are LRGs) but tends to 
over estimate the masses of the blue companions. This bias leads to an underestimation 
of merger time-scales and hence an overestimation of the fractional accretion rate for blue 
companions. 

The solid lines in Figure [3] show the calculated fractional luminosity growth of the LRGs 
assuming the dynamical friction time-scale td y n,i for the mergers. The per-subsample merger 
rates have been divided by the absolute-magnitude bin width so that the total fractional 
accretion rate is the area under (integral of) the curves. If the results in Figure [3] are naively 
interpreted as fractional mass accretion rates (they are fractional luminosity rates), the blue 
galaxies are doubly overestimated, because both the merger rate (inverse timescale) and the 
delivered mass have been over-estimated. 

Under the orbital time-scale assumption, the growth curve in Figure [3] peaks near the 
magnitude of L* galaxies. This represents fact that most of the light in the Universe, even 
near LRGs, is in L* galaxies. Under the dynamical friction assumption, the curves shift to 
more luminous, more massive, galaxies; lighter galaxies linger around the LRG for a longer 
time. 

The maximal fractional luminosity accretion rate (the sum of the integrals under the 
dashed curves in Figure [3]) is [5.6 ± 0.2] percent Gyr -1 , but this rate is unrealistically high; 
certainly pairs with large mass differences do not merge in an orbital time! The dynamical- 
friction rate is [1.7 ±0.1] percent Gyr -1 , and is also probably an overestimate because at least 
some physical pairs are not on the path to rapid merging. 

Table [1] gives the derived fractional luminosity growth from every imaging subsample 
for both the orbital time-scale assumption and the dynamical friction time-scale assumption 
as a function of the luminosity and color of the subsamples. 
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4. Discussion 

We have combined Sloan Digital Sky Survey (SDSS) spectroscopic data on 23043 lu- 
minous red galaxies (LRGs) with SDSS imaging data on enormous subsamples of fainter 
galaxies to measure cross- correlations. We have measured the projected two-dimensional 
cross-correlation functions w s i(r p ) on very small scales (0.01 < r p < 8 h^ 1 Mpc) between 
spectroscopic LRGs ("s") with luminosities [Mo.-n — 5 log 10 h] pa —22.75 mag and many sub- 
samples of imaging galaxies ("i") with luminosities [Mo.sj — 5 log 10 h] < —18 mag. The 
imaging limit is 50 times or 4.25 mag fainter than the mean LRG; the samples of companion 
galaxies cover a broad range in color and magnitude. In addition, the large volume of the 
SDSS LRG sample allows us to cut the companion galaxies into many distinct subsamples 
with different luminosities and colors but nonetheless measure the clustering as a function 
of these properties with high signal-to-noise. The principal limitation arises from the lack of 
spectroscopic information on the companion galaxies; this makes it impossible to precisely 
measure the real-space number densities for the companion subsamples. We cannot disen- 
tangle the clustering power from the number density; we only measure the product n» w S i(r p ) 
but not either or w S i(r p ) separately. 

Figures [T] and [2] show the results of these measurements of riiW S i{r p ) for red and blue 
galaxies respectively. In both figures several characteristic transition scales are visible. The 
sharp break at at r p pa 2 h^ 1 Mpc and the less-sharp transition at r p pa 0.3 h^ 1 Mpc in the 
curves can be explained in the context of the "halo occupation" picture of g alaxy clustering 



(jPeacock &: Smith! l2000t IScoccimarro et al.l l200ll ; iBerlind fc Weinberg! |2002| ) . If galaxies are 



residing within dark matter halos then the clustering of the galaxies on scales larger than 
halos is determined by the clustering of the dark matter halos that host them, plus statistics 
of the occupation of halos by galaxies. In this picture, the first transition at r p pa 2 h~ l Mpc 
locates the size of the largest halos that host LRGs — the largest halos in the Universe. At 
larger separations, at r p > 2 h^ 1 Mpc, this is the regime in which all LRG-galaxy pairs come 
from two separate halos (the "two-halo" regime). Inside this scale, at 0.3 < r p < 2 h^ 1 Mpc, 
the galaxy-LRG pairs are a mix of pairs, in some of which the companion galaxy belongs to 
the same halo as the LRG and in some of which the galaxy belongs to a separate halo. This 
"mixed" regime comes from the fact that LRGs reside in a range of halo sizes. The inner 
limit of this regime is at r p pa 0.3 h -1 Mpc, depending on the luminosity of the imaging galaxy 
subsample i in question. This inner scale is close to the virial size of the smallest halo that 
can host an LRG (plus the virial size of the smallest halo that can host a companion galaxy 
from the imaging subsample i). At smaller scales, at r p < 0.3 h~ l Mpc, all the galaxies belong 
to the same halo as the LRG halo (the "one-halo" regime). Here the clustering represents 
the mean radial profile of the halos mixed with details of galaxy evolution. Figure [1] shows 
an increase in the clustering of blue galaxies at r p pa 50 h~ l kpc toward the central regions of 
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the halo. This could be caused by a boost of star formation in these galaxies, which makes 
them more luminous and places them in higher luminosity bins in our calculation. 

Finally, both Figures show a sharp drop in the clustering power on scales r p < 30 h -1 kpc. 



This could be due to failure of the object detection software of the SDSS (e.g.. iMasjedi et al. 



20061 ). or it could be a real effect from disintegration of galaxies by dynamical friction or 



other tidal stripping expected in some galaxy evolution models. 

We integrate the de-projected, three-dimensional cross-correlation functions rii£ S i(r) for 
each imaging subsample i on very small scales to calculate the average number of galaxies 
that are in dynamical pairs with each LRG. We use two different time-scales for merger 
events to calculate merger rates. The first is the orbital time t m ut, equivalent to assuming 
that all galaxies merge in one orbit. This is the shortest time imaginable to merge, so it 
provides a strict upper limit on the merger rate. The second time-scale is the dynamical 
friction time-scale idyn,i for which we approximate the merger time with a linear function of 
the mass ratio (m s ) / (mi) of the mean galaxies from the two samples. This is equivalent to 
assuming that equal-mass (LRG-LRG) mergers take the one-orbit time, but pairs of galaxies 
with more different masses take longer times to merge. 

We use the two time-scales to calculate both a strict upper limit to the merger rate and 
a more realistic rate, although even the dynamical friction time-scale calculation involves 
assuming that essentially all close pairs merge. For both time-scales we have measured the 
fraction of LRG luminosity that is added to the LRG through mergers of galaxies from each 
imaging subsample % per Gyr. The fractional luminosity growth for LRGs through mergers 
as a function of the color and magnitude of the merging companions is shown in Figure [3j 

Most of the luminosity brought into LRGs by merging is brought by red companions or 
"dry mergers," and most of it is brought by galaxies near (or above) L*. The contribution to 
growth decreases with decreasing luminosity at the faint end; the curves essentially to zero 
by [Afo.3j — 5 log 10 h] = —18 mag. Calculation of the total amount of luminosity brought in 
by merger activities does not require consideration of fainter companion galaxies. 

Integration of Figure [3] over companion absolute magnitude yields the total fractional 
amount of LRG luminosity growth from mergers. Under the maximal one-orbit merger time 
assumption, we find that the total fractional growth rate is strictly smaller than [5.6 ± 
0.2] percent Gyr -1 . Under the more realistic dynamical friction assumption, we find that 
the total fractional growth rate is [1.7 ± 0.1] percent Gyr -1 where only about one-tenth of 
that is through "wet mergers" (blue companions) and the rest is through dry mergers (red 
companions) . 

Our results are not consistent, at face value, with most morphological measurements of 
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the merger rate — measures that involve identification of merging galaxies by their appearances 
most of w hich find rates on the order of ten percent per Gyr or of order un i ty over a Hubble 
time (e.g., lAbraham et al.lll996t ; IConselice et all 120031 : Ivan Dokkuml 120051 : lLotz et all 12006 



though see also lDe Propris et alll2007l ). Our method for inferring the merger rate suffers — as 
all these other investigators' methods do — from uncertainties in merger timescales. However, 
we avoid all issues related to morphological selection of merging systems or merger remnants, 
which tend to introduce subjectivity, and we avoid uncertainties related to line-of-sight pro- 
jections because we work with the true three-space correlation function. Our method is 
therefore much more precise than other methods. Furthermore, because merging cannot 
happen on timescales shorter than a dynamical time, our upper limit is extremely robust 
(and not made uncertain by projection effects). 

Our merger rate estimate can be reconciled with other estimates if we assume that either 
(1 ) the merger rate is an extremely strong function of the primary galaxy mass (since we only 
investigate the rate for the most massive galaxies in the Universe), or (2) merging produces 
observable distortions to galaxy morphologies (e.g., tidal tails) that last for many dynamical 
times, or (3) significant morphological signs of merging can be raised by very frequent, very 
minor mergers, which don't contribute much to the build-up of mass. Our r esults are more 



consi s tent with measures of the merger rate based on counts of close pairs ( Carlberg et al. 



1994; 


Patton et al. 


1997 




van Dokkum et al. 


1999; < 


Carlbere; et al. 


2000: 


Patton et al. 


2000; 


Lin et al.ll2004: 


Masjedi et al. 


2006; 


Bell et al 


2006< 


»), but even there, our results are on the 



low side. 



Murali et al. 


2002; 


Mailer et al. 


2006 





20071 ). but we caution that no predictions have been made for exactly what we have observed, 
and that galaxy-galaxy merging occurs at length and dynamical scales where cosmological 
simulations are not completely reliable. 

There are three respects in which the luminosity growth shown in the solid lines in 
Figure [3] — the more "realistic" estimates — are nonetheless over-estimates or upper limits on 
the true fractional mass growth for LRGs. First, we are assuming that the vast major- 
ity of pairs do merge as quickly as dynamical friction allows. This is not true for close 
pairs in high velocity-dispersion environments. In ad dition, even when the pairs are bound 
the Chandrasekhar formula may be an overestimate (jBoylan-Kolchin et all 120071 ). Second, 
the blue galaxies have both their masses over-estimated and their dynamical friction times 
under-estimated with a constant mass-to-light ratio assumption, so the blue galaxies do not 
contribute as much mass as Figure [3] implies. Third, we are assuming that all the stars in 
the companion galaxies will end up in the central LRGs. Recent work has suggested that 
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this is not the case and in fact up to 50 percent of the stars in the companions could be 
stripped off the companion before the merger is complete and contribute to the intra-cluster 
light instead of the luminosity of the LRG (ILin et all 120041 ) . This suggests that even our 
dynamical friction assumption could still be an upper limit on the growth of the LRGs. 

These results are consistent with recent results on the evolution of the luminosity func 



Blanton 


2006; 


Wake et al. 


2006; 


Faber et al. 


2007; 


Brown et al. 


2007) 



at face value and assume that the growth happens at a non-evolving rate, we expect the 
LRGs to grow by about « 10 percent between redshift z — 1 and z = 0.1 (a period of 
« 6 Gyr) . 
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Fractional luminosity growth of LRGs from mergers 
[Mo.3j — 5 log 10 h] Red max Blue max Red DF Blue DF 
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Table 1: Measurements of the fractional growth of spectroscopic LRGs over a Gyr through 
merger events with different imaging subsamples of companion galaxies with different mag- 
nitude ranges and colors. The second and third ( "max" ) columns present the measurements 
under the maximal assumption that all pairs merge in an orbital time. The fourth and fifth 
("DF") columns present the measurements under the more realistic dynamical-friction as- 
sumption. Note that the "max" and "DF" columns are given in units that differ by a factor 
of 10. 
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Fig. 1. — Projected two-dimensional cross-correlation functions UiW s 



of spectroscopic 



LRGs s with red companion imaging galaxy subsamples % with different luminosity ranges, 
weighted by the number density of each imaging subsample, and scaled by r p for better 
illustration. The spectroscopic LRG subsample s has absolute magnitudes —23.2 < [Mo.s g — 
5 log 10 h\ < —21.2 mag and has been trimmed to redshifts 0.16 < z < 0.30. The uncertainties 
are estimated by jackknife, with jackknife trials dropping contiguous sky regions (see text). 
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Fig. 2. — Same as Figured] but for blue companion imaging galaxy subsaples i. 
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Fig. 3. — The mean fractional luminosity growth of LRGs per Gyr per unit absolute mag- 
nitude of the companion, derived from the data in Table [TJ The dashed thick red and thin 
blue lines show this quantity for the red and blue companions respectively, made under the 
maximal assumption that all companion galaxies merge into the LRG in one orbital time. 
The solid thick red and thin blue lines show the same thing but made under the more real- 
istic assumption that companion galaxies merge in a time governed by dynamical friction. 
The total fractional growth rate is the integr al (area) under the c urves. The green hatched 
region shows the result of our previous work (IMasjedi et al.l 120061 ) . 



